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ABSTRACT 



Context, FK Comae is a rapidly rotating magnetically active star, the light curve of which is modulated by cool spots on its surface. 
It was the first star where the "flip-flop" phenomenon was discovered. Since then, flip-flops in the spot activity have been reported 
in many other stars. Follow-up studies with increasing length have shown, however, that the phenomenon is more complex than was 
thought right after its discovery. 

Aims. Therefore, it is of interest to perform a more thorough study of the evolution of the spot activity in FK Com. In this study, we 
analyse 15 years of photometric observations with two different time series analysis methods, with a special emphasis on detecting 
flip-flop type events from the data. 

Methods. We apply the continuous period search and carrier fit methods on long-term standard Johnson-Cousins V-observations from 
the years 1995-2010. The observations were carried out with two automated photometric telescopes, Phoenix-10 and Amadeus T7 
located in Arizona. 

Results. We identify complex phase behaviour in 6 of the 15 analysed data segments. We identify five flip-flop events and two cases 
of phase jumps, where the phase shift is Acf> < 0.4. In addition we see two mergers of spot regions and two cases where the apparent 
phase shifts are caused by spot regions drifting with respect to each other. Furthermore we detect variations in the rotation period 
corresponding to a differential rotation coefficient of \k\ > 0.031. 

Conclusions. The flip-flop cannot be interpreted as a single phenomenon, where the main activity jumps from one active longitude 
to another. In some of our cases the phase shifts can be explained by differential rotation: Two spot regions move with different 
angular velocity and even pass each other. Comparison between the methods show that the carrier fit utility is better in retrieving slow 
evolution especially from a low amplitude light curve, while the continuous period search is more sensitive in case of rapid changes. 

Key words. stars:late-type - activity - starspots - dynamo - HD 1 17555 
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1. Introduction 

FK Comae Berenices (HD 117555; hereafter FK Com) is the 
prototype of a class of single chromospherically active rapidly 
rotating G-K gi ants. Only a few stars fu lfil the definition of the 
FK Com class (Bopp & Rucinskil ll98ll) . These stars may rep- 
resent an interm ediate state of coalesced W UMa binaries (e.g. 
Bop p & Stencell 19811) in the process of magnetic braking, which 
would explain why they are so rare. 

FK Com itself is an extremely active late-type star and has 
been extensively studied with ground-based optical photome- 
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* Based on data obtained with the Amadeus T7 Automatic 
Photoelectric Telescope (APT) at Fairborn Observatory, jointly op- 
erated by the University of Vienna and AIP, the Phoenix-10 
APT at Mt. Hopkins, Arizona, and the Nordic Optical Telescope, 
Observatorio Roque de los Muchachos, La Palma, Canary Islands. 
The photometric observations are available in electronic form at the 
CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http://cdsarc.u-strasbg.fr/vizbin/qcat ?J/A+A/yyy/Axxx 



try and spectroscopy, radio observations, a s well as s atellite- 
based UV- and X-ray observations (e.g . IJetsu et al.1 1994a 



llet; 



2009a; Hughes & McLea 



l Olah etalJ 120061; iKorhonen et al.l . 

1987: Bopp & Stencell 1 198 It lAvres et all 120061 iDrakeetaLl 
2008). 



The photomet r ic rotation period of FK Com is P 
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(IChugainovlll966t IJetsu et al.lll993|). It has been proposed tha t 
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its spectra l class is between G5 m dKorhonen et al.ll 1 999L l2007h 
and G4 in (Strassmeier 2009). Korhonen et al. (2000) concluded 
that v sin i = 159 km s _1 gave the best fit for the spectral data. 
Analysing photometry spa nning roughly over 25 years 
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(1966-19901. 1 Jetsu et all (1 19931) reported a switch of activity be- 
tween two longitudes separated by approximately 180°, and la- 
belled this effect the "flip-flop". The activity was observed to 
jump from one active longitude to the other three times during 
the period of analysis. Furthermore, the active longitude system 
was reported to be rotating with the photometric rotation period 
of Pphot = 2. d 4002466 + 0. d 0000056 throughout the whole span 
of the data. 

The flip-flop phenomenon has since then been sug- 
gested to occur in a number of stars (e.g. Jetsu 1996; 
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iBerdv ugina & Tuo minenll998tlLehtinen et al.l2011h . While fur- 
ther results of the stable active longitu de system p r oducin g flip- 
flops on FK Com were published by [Jetsu et all (Il994bh . evi- 
dence for the phenomenon being more com plex started build- 
ing up, e.g. by the analysis of photometry by Strassmeie r et alj 
d!997al) . In this study, a gradual drift of the spots from one ac- 
tive longitude to the other was detected during 1993-1995, in 
contrast to the abrupt changes reported earlier, where spot mi- 
gration over phase was not related to the phenomenon. With im- 
proved photometric data with denser timing piling up, the picture 
of a steady active longitude system rotating with one single pe- 
riod w as abandoned. For instance in the study of lKorhonen et alj 
(I2002h . where one-dimensional photometric inversions of FK 
Com were presented, the active longitude system exhibited three 
different periods: in the beginning of the observations, the sys- 
tem was rotating with the photometric rotation period of the star, 
slowed down d uring 1994-1997, and sped up to super-rotation 
for 1998-2004 dKorhonen et alj|2004l) . The first apparent semi- 
regularity of the flip-flops also became under doubt; the time 
between the flip-flops could range from a year to several years. 
Cycles of 5.2 and 5.8 years in the migration o f the two active lon- 
gitudes were reported by lOlah et all (2006). Furthermore they 
made a distinction between phase jumps and flip-flops. In the 
case of phase jumps, a new active region appears on the same 
hemisphere of the star as the old active region, resulting in a 
phase shift of less than 0.5. During a flip-flop the active longi- 
tude changes ~ 180°. 

The behavior of the active longitude system and the re- 
lated flip-flops were spectroscopically confirmed by analysing 
the uniquely long series of Doppler imag es from the years 1993- 
2008 by Korhonen and collaborators dKo rhonen e t alJ I2000L 
120071: lAvres et alj 120061: iKorhonen et al l l2009allbl) . "This work 
was based mainly on the observations with the SOFIN high- 
resolution spectrograph at the Nordic Optical Telescope (La 
Pal ma, Spain) 

iJetsu et all d!994bl) detected variations in the photometric pe- 
riod, which could be a signa ture of differential rotation in FK 
Com. IKorhonen et alj ((2007) estimated the differential rotation 
by combining Doppler images and period analysis of photome- 
try. They reported a rotation law of 

Q ~ (151.307d + 0.09/d)-(1.787d±0.12 o /d)sin 2 ^ (1) 

where i/r is the stellar latitude. The estimated relative differential 
rotation coefficient was thus k ss 0.012. It should, however, be 
emphasized that estimating differential rotation using Doppler 
imaging is challenging, especially because of artifacts and errors 
in the spot latitudes. 

Observations of FK Com obtained with the Far Ultraviolet 
Spectroscopic Explorer (FUSE) revealed complex profiles sug- 
gesting that the transition region and the corona are highly 
structured , dominated by dy namic processes, and could be very 
extended dAvres et alJ l2006). The X-ray observations obtained 
with XMM-Newton imply that the corona of FK Com is dom- 
inated by large magnetic structures similar to the intercon- 
necting loops in so lar active regions, but significantly hotter 
dGondoin et al.l2()02[). Similarly, us ing data from the Chandra X- 
ray Observatory, iDrake et all d2008l) found indications of mag- 
netic loops in the corona of FK Com. Their observations sug- 
gested that the observed X-ray emission originates from plasma 
residing predominantly in extended structures centred at a phase 
halfway between two spot regions, and that the coronal struc- 
tures revealed by the Chandra observations correspond to mag- 
netic loops joining these two spot regions. This would support a 
model where the two regions have opposite magnetic polarities. 



Table 1. Summary of the observations and labelling of each seg- 
ment in the text. 



Segment 


^rnin 




Telescope 




HJD-2400000 


HJD-2400000 




SEG1 


49876.7555 


49908.6573 


PhlO 


SEG2 


50085.0551 


50265.6741 


PhlO 


SEG3 


50412.0383 


50636.6781 


PhlO 


SEG4 


50778.0423 


50997.7057 


PhlO &T7 


SEG5 


51144.0394 


51362.7280 


PhlO & T7 


SEG6 


51508.0427 


51731.6859 


PhlO & T7 


SEG7 


51873.0421 


52089.7274 


PhlO & T7 


SEG8 


52242.0328 


52461.6984 


PhlO & T7 


SEG9 


52613.0093 


52828.6937 


PhlO & T7 


SEG10 


52972.0288 


53194.6817 


PhlO & T7 


SEG11 


53343.0462 


53565.6808 


T7 


SEG12 


53709.0464 


53922.6978 


T7 


SEG13 


54075.0453 


54285.6861 


T7 


SEG14 


54440.0490 


54643.7298 


T7 


SEG15 


54807.0399 


55004.7444 


T7 


SEG16 


55172.0443 


55297.8314 


T7 



We want to investigate the spot activity in more detail and 
combining different methods. Our main aim is to identify and 
study the flip-flops of FK Com. We apply two novel time se- 
ries analysis methods on long-term p hotometric observa tions; 
the continuous period search method (Lehtinen e t al.l201 U here- 
after CPS) and the carrier fit utility dPelt et al.ll201 ll hereafter 
CF). In particular, we are interested in the nature of the flip-flops. 
We want to isolate different types of proposed effects: i) abrupt 
jumps from one hemisphere to another representing the flip-flop 
phenomenon, ii) events better described as phase jumps (with a 
phase separation significantly less than 180°), and iii) gradual 
phase drifts of the active longitudes. We also aim at establishing 
whether, when properly classified and separated, there are any 
regularities related to the phenomena. 

Our secondary aim is to compare the results from our two 
analysis methods. We anticipate that there is a great advantage of 
combining the two time series analysis methods because of their 
different solutions for modelling the data. It is also important to 
make a comparative study with real data. 

2. Observations 

The photometric standard Johnson-Cousins V observations 
were collect ed with the 0.75m Vie nna University/AIP APT 
"Amadeus" dStrassmeier et alj|1997bl T7 in TablsQ}, located at 
Fairborn Observatory and the Phoenix- 10 APT at Mt. Hopkins, 
Arizona (PhlO). Data points with errors greater than 0.02 were 
rejected and the differential magnitudes were transfered into ap- 
parent m agnitudes u s ing the the comparison star HD 117567 
just as in IJetsu et all dl994al) . For a more detailed description 
of the observat i ons o f FK Comae, we refer to the paper by 
IKorhonen etail d2001l) . 

Most of the data has been included in earlier papers 
( Kor honen et all [2002h lOlah et al.ll2006l 120091 IKorhonen et all 
2009a, and references therein). The data is summarised in 
Tableland published electronically in the CDS. Note that each 
observing season forms a segment. 

3. Analysis methods 

Both the CPS and CF analysis methods are based on an approach 
of continuous curve fitting. A main difference is that the CPS 
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Fig. 1. All photometric V-data for FK Com. The Phoenix- 10 data is marked with black dots and the Amadeus data is marked with 
grey plusses. 



method allows the period to vary, while the CF uses a constant 
carrier period. Both approaches can be argued for. On one hand 
e.g. differential rotation will cause the photometric period to de- 
pend on the spot latitude or the anchoring depth of the spot struc- 
ture. On the other hand, noise in the data and insufficient phase 
coverage will contribute to spu rious period variations in the CPS 
method dLehtinen et al.ll201ll) . The CF method is more stable 
against such errors. 

3.1. CF method 

The CF method is based on the simple idea of decomposing the 
observed stellar light curves into two components: 1) a rapidly 
changing carrier modulation tracing the regular part of the sig- 
nal, for instance rotation of a spotted star, and 2) a slowly chang- 
ing modulation, such as evolution of the cool spots on the stellar 
surface. Such a situation can be described with the following 
mode 

K 

y ct (t) = fl (f) + Yj (a k (t) cos(2**/ f) + b k (t) %m{2nkf f)), (2) 

where ao(f) is the time-dependent mean level of the signal, K is 
the total number of harmonics included in the model, describ- 
ing the overtones of the basic carrier frequency, while a^(f) and 
bk(t) are the low-frequency signal components. The carrier fre- 
quency fo can be either known a priori, or determined using the 
CF utility as the first step of the analysis. In this paper, we take 
the previous determinations of the photometric rotation period as 
the first guess of the carrier period Pq = l/fo - 2 d .40. The next 
step in the analysis is to formulate a suitable model for the mod- 
ulating curves. In lPelt et alJ JIoTlh we introduced two classes of 
models based on either trigonometric or spline approximation. 
In this paper, models based on the trigonometric approximation 
are used. 

The trigonometric approximation model for the modulat- 
ing signals is built in the following way. Let the time interval 
[tmin, Wx] be the full span of our input data. Then we can in- 
troduce a certain period D = C x (f max - f m ; n ) for which the 
coverage factor C is larger than unity (typically C = 1.1 - 1.5). 



Using the corresponding frequency, fo = l/D, we can now build 
a trigonometric (truncated) series of the type: 

L 

a(t) = eg + V (c a t cos(2ntlf D ) + s", sm{2ntlf D j), (3) 
l=i 

and 

L 

b(t) = 4 + J] (cf cos{2ntlf D ) + sf sm(2ntlf D )), (4) 

where L is the total number of harmonics used in the modulator 
model. According to our definition of a slow process, the pe- 
riod D must be significantly longer than the carrier period Pq. 
With the data segment lengths in the regime of 100-200 days, 
this condition is well satisfied with the chosen coverage fac- 
tor. Next, proper expansion coefficient estimates are computed 
for every term in the series for the fixed carrier frequency fo 
and the "data frequency" fy; this is a standard linear estimation 
procedure and can be implemented using stan dard mathematica l 
(statistical) packages, as described in detail bv lPelt et al.l d201 lb . 
If all coefficients {a k , b^) consists of the same number of har- 
monics L and we approximate separate cycles by a /f-harmonic 
model, then the overall count of linear parameters to be fitted is 
Af = (2xL+l)*(2xA'+l). The actual choice of the represen- 
tative parameters K and L depends on the particular object we 
are working with. The number of tones, K, depends on the com- 
plexity of the phase curves. The choice of L is constrained by 
the longest gaps in the time series. In this study, we adopt K — 2 
and L = 3, resulting in the total number of free parameters to be 
fitted N = 35. 

We visualise our results in the following way. First we cal- 
culate a continuous curve least-squares estimate y c {(t), from the 
randomly spaced and gapped data set. This approximation is 
continuous and does not contain gaps, and therefore allows us 
to get a smooth picture of the long-term behaviour. Next we di- 
vide this continuous curve into strips with a length of the carrier 
period Pq = jr. We then normalize each strip so that the approx- 
imating values span the standard range of [-1, 1]. After normal- 
isation, we stack the strips along the time axis. To enhance the 
obtained plot, we extend every strip somewhat along phases, so 
that the actual display is wider (along phases) than a single pe- 
riod. The normalisation is a relevant part of our procedure be- 
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cause it helps to grasp the phase information we are interested in 
(trends, drifts, flip-flops etc). This method of visualisation allows 
to verify that the model fits into the data and not into the gaps. If 
the number of nodes or harmonics L for the modulation model 
curves is chosen properly, then the phase plots do not reveal any 
underlying timing structure. 

The data from both telescopes were merged for the CF- 
analysis. In principle merging V-magnitudes from two different 
sources may pose problems, even though the reduction is care- 
fully done. In this case however, we concluded that the data was 
uniform enough for the CF-analysis. The CF analysis was ap- 
plied on segments SEG2-16; in SEG1, there were not enough 
points to carry out the analysis. 

The goodness of the fit was estimated with the coefficient of 
determination 

i i 

where y is the mean of the data. Our goal was to have R 2 > 0.9. In 
five out of all the fifteen segments R 2 was too small. Therefore, 
for SEG4, SEG5, SEG9, SEG10 and SEG15, we refined our 
analysis by excluding the 3cr outliers after an initial CF fit, and 
making a new fit. It is probably no coincidence, that four of these 
segments consisted of data from the two different telescopes. 
However, tests showed that excluding the outliers had no other 
significant influence on the result than that of improved good- 
ness. 

3.2. CPS method 

The CPS method was originally developed for the analysis of 
photometric observations of late-type stars. In order to under- 
stand the spot activity of these stars, one has to take into ac- 
count both short time scale (days) and long-term (years) changes 
( Leht inen et al.l l201lh . The method is base d on the Three Stage 
period Analysis (TSPA. lJetsu & Peljl 19991) . The model 

A' 

y cps (ti,P) = ao + 2^ \- a k cos (Znkftd + b k sin (2n*/f,-)], (6) 

k=\ 

is used to fit each set of the data. Here /? consists of the param- 
eters (a,, hi, /), which are determined through a non-linear least 
squares optimisation. The three major improvements compared 
with the TSPA are: 

1 . We analyse the data with a sliding window of length Ar max 
in order to increase the time resolution. 

2. We test models of different orders and choose the final order 
K using a Bayesian information criterion . 

3. We derive the time scale Tc of significant changes in the light 
curve. 

With the CPS-analysis one can thus derive a continuous se- 
ries of estimates for the mean magnitude (M), total light curve 
amplitude (A), photometric period (P) and epochs of the light 
curve minima (f m ; n ). The temporal changes of the mean and am- 
plitude are useful for studying variations in the level of spot ac- 
tivity, since the mean magnitude will be sensitive to the spot cov- 
erage and the amplitude is a measure of the non-axisymmetry of 
the spot configuration. Variations in the photometric period can 
be caused by differential rotation, or alternative ly dynamo waves 
dKrause & Radlerill980t iTuominen et al.ll2002l) . The time scale 
of significant change Tc can be used to estimate the stability of 



the light curve. It may also be relat ed to the convective turnover 
time r c dLehtinen et al.ll201 lLl2012l) . 

The CPS-analysis was applied to the Phoenix 10 and 
Amadeus data separately. This method uses less data points in 
each fit than the CF-method, which makes it more vulnerable to 
errors in the data. Furthermore, with K = 2 the number of free 
parameters is six, i.e. much less than in the CF-method. Thus, 
the need to maximize the number of data points is less impor- 
tant than the homogeneity of the data for the CPS analysis. The 
maximum length of the moving window was Ar max = 24d. This 
length was chosen because the rotation period being ~ 2?4, it 
will give an optimal phase coverage in the case of evenly spaced 
observations. Normally Ar max would als o define the div i sion o f 
the data into segments, as described in Lehtinen et al] d201 lb . 
However, for consistency reasons in this paper the segments are 
identical to the observing seasons. The maximum order of the fit, 
i.e. the largest tested /f-value for the model (Eq.|6]l, was Kn m = 2. 

3.3. Kuipertest 

Non-parametrical time-series an alysis method s, such as the 
Kuiper test for phase distributions (Kuiper 1960), can be utilised 
to identify active longitudes from the epochs of light curve 
minima. We used t he unweighted Kuiper-test as formulated by 
iJetsu & Pelt! {T996). The Kuiper periodogram is calculated for a 
set of epochs of photometric minima t m \ n . The most significant 
periods are tested against the null hypothesis of a random phase 
distributio n. Exampl e s of application of t his method can be 
found in Hetsul d 1 9961) . iLehtinen et all d20 1 ll) and lLehtinen et all 
d20Tl . 

We computed the Kuiper statistic periodogram for period- 
icities of 2.2 - 2.6 days using the epochs of photometric min- 
ima f m ; n derived both with the CF and CPS methods. A total 
of 1637 primary and secondary minima where retrieved from 
the CF-analysis. From the CPS-analysis we used separately the 
136 independent minima and all 1837 minima. In addition to the 
minima, we also tested 1640 maxima from the CF-analysis. 

4. Results 

The results from the CF analysis are shown in the fifteen panels 
of Fig. [2] wrapped with the carrier period Pq = 2 < ?40 and as- 
signing the phase <p — to the first time point of each segment. 
The phase information of the segments is therefore not compara- 
ble. The dark-bright pattern is replicated as the phase-axis is ex- 
tended around [0, 1] to help the visualisation. The dark colours 
represent higher magnitude, i.e. lower temperatures. The results 
of the fitting procedures are summarised in Table [2] where R 2 
and R\ are the coefficients of determination (Eq. [5]) before and 
after the removal of the 3<x outliers, respectively. AN is the num- 
ber of removed data points. As an example, the fits for two seg- 
ments are shown in Fig. [3] 

The full CPS results of segments SEG4 and SEG5 are shown 
in Figs.[4]and[5] The values of the a priori period estimate Pq, the 
median of all reliable periods P me d, the limiting modelling order 
Ajj m and the maximum length of the dataset Ar max are given at 
the top of the plots. The panels show: 

(a) standard deviation of residuals <x £ (t), t being the mean 
epoch of each data set; 

(b) modelling order K(t) (squares, units on the left y-axis); and 
number of observations per dataset n (crosses, units on the 
right y-axis); 

(c) mean V-magnitude M(t); 
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Segment 02 Segment 03 Segment 04 




50 100 150 50 100 150 50 100 150 
Segment 05 Segment 06 Segment 07 Segment 08 




50 100 150 50 100 150 50 100 150 50 100 150 
Segment 09 Segment 10 Segment 1 1 Segment 12 




50 100 150 50 100 150 50 100 150 50 100 150 




Time (days) Time (days) Time (days) Time (days) 



Fig. 2. CF analysis results for each segment with the carrier period Po = 2 d .40. Each panel shows the computed time-dependent 
phase diagram, i.e. the light curve amplitude profile over phase (y-axis) plotted as function of time (x-axis). 




50800 50850 50900 50950 51150 51200 51250 51300 51350 

HJD-2400000.0 HJD-2400000.0 



Fig. 3. The CF fits for segments SEG4 and SEG5. The residuals at the bottom of the plot are shifted with 8.28 (SEG4) and 8.23 
(SEG5). 



(d) time scale of change 7c(t); 

(e) amplitude A(r); 

(f) period P(r); 

(g) primary (squares) and secondary (triangles) minimum 
phases m i n ,i(T) and 4>min.i{ T )- These phases are calculated 
using the median period P me d of the segment; 

(h) M(t) versus P(r); 



(i) A(t) versus P(t); 
(j) M(t) versus A(t). 

In the subplots (a), (c) and (e)-(g), the reliable parame- 
ter estimates are indicated by filled symbols and unreliable 
ones by open symbol s. The reliability is tested as described by 
iLehtinen et al.l d201 ll) . In subplot (d), the upward pointing ar- 
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FK Comae, SEG = 4 
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Fig. 4. The CPS analysis of segment SEG4 of the Amadeus photometry of FK Com. Descriptions of the subplots are given in Sect, 
i 



Table 2. Summary of CF fit results for each segment. 



Segment 


R\ [%] 


AN 


K [%] 


SEG1 








SEG2 


95.3 






SEG3 


96.4 






SEG4 


62.9 


4 


83.0 


SEG5 


83.6 


6 


86.7 


SEG6 


95.0 






SEG7 


97.7 






SEG8 


93.0 






SEG9 


88.0 


6 


90.7 


SEG10 


78.5 


6 


89.9 


SEG 11 


95.7 






SEG 12 


97.1 






SEG 13 


98.0 






SEG14 


93.1 






SEG 15 


89.1 


2 


91.0 


SEG 16 


98.5 







rows signify that no significant change was detected before the 
end of the segment, i.e. the model was valid from this point 
throughtout the segment. 

In the correlation plots (h)-(j), the error bars have been 
drawn only for the independent parameter estimates. The linear 
Pearson correlation coefficients ro for the independent datasets, 
as well as an estimate of the probabilities P(\r\ > ro), are given. 

All independent CPS period estimates are shown in Fig. [6] 
The mean period and its standard deviation were P w + AP W w 



2^3975 + $0123. There is no doubt, that the photometric period 
is varying. The probable reason is differential rotation. The vari- 
able period poses problems on what to use as a standard period 
for plots. We use the "old" ephemeris 

HJD min = 2439252.895 + 24002466£ (7) 

derived bv lJetsu et al.l (1 19931) whenever we are combining results 
from the CF, CPS or Doppler imaging analysis. In the rest of the 
plots, we use either the carrier period (Fig. |2}, the CPS median 
period for each segments (Figs.|4]and|5]l or the best active longi- 
tude period (Fig.|7]i. 

The independent CPS estimates of the mean magnitudes and 
light curve amplitudes are plotted together with the correspond- 
ing results from the CF analysis in Fig [6] In Fig. [8] we show 
the CF and CPS photometric minima for some interesting seg- 
ments, together with results from Doppler imaging when avail- 
able. In the latter cases we plot latitudinally averaged slices of 
Doppler images calcu l ated by H . Korhonen (lAvres et al.l l2006: 
iKorhonen et al.l 120071 l2009allbl) . These slices were calculated 
from the original rectangular DI maps by weighting each sur- 
face element with its size. For segment SEG5 we also plot all 
reliable CPS light-curves in Fig. [9] 

We can see that there is no discrepancy between the CF and 
CPS methods, since all differences can be explained by the dif- 
ferent approaches. The results from the CF method look like a 
smooth fit to the more noisy CPS results. 

The mean of Tc(t) is 32 d and the minimum 7 d. This is both 
surprising and worrying. In previous analysis of main sequence 
stars Tc has been of the same order as the convective turnover 
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Fig. 5. The CPS analysis of segment SEG5 of the Amadeus photometry of FK Com. 



time r r (iLehtinen et alJl201ll [2Q12h . It is hard to determine r r 
for FK Com since it is an exceptional giant star with uncertain 
parameters. However, with an estimated radius of R ~ IORq, 
one would expect t c to be larger than for main sequence stars. 
Thus one would expect also Tc(t) to be much larger. The wor- 
rying part is that about 30 % of the Tc(j) are smaller than the 
Ar max = 24 d used as the time window for the CPS-analysis. 
Tests with reducing Ar max showed, however, that there was no 
significant change in the result, except that the reduced num- 
ber of points in each CPS-set increased the errors. A low Tc(t) 
indicates that rapid changes occur in the spot configuration. It 
should be emphasised that Tc(j) of course depends on K: The 
higher the order of the model, the more sensitive it is to changes. 
Furthermore, the analysis is based on statistics, which explains 
why some low amplitude models may succeed in fulfilling the 
whole data, while subsequent models may be quickly rejected 
(see e.g. the three "arrows" in the upper left corner of panel 
(d) in Fig. [5). Deta ils on the estimation of Tc(j) are found in 
ILehtinen etallfcO lib . 

4. 1 . Active longitudes 

The overall result is that in about half of the segments only one 
active longitude can be detected. Two active longitudes usually 
seem to be present before and after longitudinal shifts in the ac- 
tivity, but there are also cases, e.g. SEG16, of two simultaneous 
active longitudes without indications of shifts. In some panels of 
Fig. [2] the pattern runs slanted upwards, in a couple it remains 
nearly horizontal, in a few panels it is slanted downwards, or 
in some segments disrupted in some way or another. The rising 



(e.g. segments SEG2 and SEG3) or falling (e.g. SEG6, SEG9 
and SEG 13) trends are indications of the carrier period not be- 
ing optimal for the segment, i.e. rising trends could be corrected 
with increasing the carrier period somewhat, and correspond- 
ingly falling trends by decreasing it. As both types of trends are 
present in the segments, however, a global CF analysis of all 
the segments together would, in any case, give a carrier period 
very close to the one already adopted. The trends are relatively 
short-lived, as they are visible only in two, maximally three con- 
secutive segments, i.e. last roughly one year. A drift longer by 
an orde r of magnitude has been reported in the RS CVn binary 
II Peg dLindborg et alj 1201 U lHackman et alj 1201 lb during the 
years 1994-2001. It was interpreted as a possible azimuthal dy- 
namo wave arising from the properties of the non-axisymmetric 
dynamo solution. During these years, the spot-generating struc- 
ture was observed to rotate with a shorter, but constant period, 
i.e. forming a more or less rigidly rotating entity. 

Therefore, the short-lived trends are likely to be related to ro- 
tational non-uniformities either on the surface or at larger depth 
that the spots may be anchored to. This can also be seen in the 
varying photometric period retrieved with the CPS method (Fig. 
[5). With the seasonal gaps in the data, these trends make it hard 
to visually follow any long-lived active longitudes in FK Coma 
(see e.g. Fig. [Tot. 

However, the Kuiper periodogram analysis of the indepen- 
dent photometric minima from the CPS-analysis gave the best 
period P. di * 2^401 151 + O?000092 with a significance level of 
Q « 5.3 • 10~ n . The phases of the CPS minima, folded with 
this period, are displayed in Fig. [7] We also applied the Kuiper 
test on all photometric minima retrieved by the CPS method, 
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as well as all minima and maxima from the CF-analysis. For 
these tests extremely low Q-values were derived, but since the 
measurements cannot be seen as independent, this has no phys- 
ical relevance. All the CPS-minima yielded the period P a \ « 
2?401 173 + $000015 and the CF-minima gave the result P a * 
2?40 11668 + $0000091. However, the most significant period 
for the CF-maxima was P max * 2^405497 + $000042. 

Thus, the analysis of the CF-minima gave practically the 
same P d \ as for the CPS-minima. We note that P d \ is slightly 
longer than the mean photometric period P w . The value of 
P a i should describe the possible period of a magnetic structure 
within the star, while P w may reflect the rotation at the surface. 
This would imply, that there is a magnetic structure rotating 
slightly slower than the surface of the star. The best period for 
the CF-maxima P mdx was again slightly longer than P d \. Thus, 
this periodicity may describe something else than the mere ab- 
sence of spots. 

4.2. Flip-flop like events 

Switching active longitudes were detected in six of the 15 anal- 
ysed segments, namely SEG4, SEG5, SEG8, SEG11, SEG12 
and SEG15. These events are summarised in Tableland plotted 
in Fig. |8]together with longitudinal slices of Doppler images. In 
many cases these events were not "proper" flip-flops. Either the 
phase jump was considerably less than 0.5, or the event involved 
more spot evolution, than just a shift from one active longitude to 
another one. In order to make a distinction between these events 
we use the definitions for a flip-flop suggested by Kajatkari et al. 
(2012, A&A, in prep.): 

- the region of main activity shifts about 180 degrees from the 
old active longitude and then stays on the new active longi- 
tude. 

or 

- the primary and secondary minima are first separated by 
about 180 degrees, after which the secondary minimum 
evolves into a long-lived primary minimum, and vice versa. 

In SEG4 (HJD = 2450778.043), two phase changes can be 
detected from the CF analysis (upper left panel of[8]l. During the 
first phase disruption, the primary minimum segregates into a 
secondary minimum roughly at HJD a 2450870. The two min- 
ima have a phase separation of A0 » 0.4, which is close enough 
to 0.5 to call the event a flip-flop. At around HJD a 2450880, the 
secondary minimum has become the main minimum. The CPS 
analysis shows, that the two parallel minima persist for approx- 
imately 35 days. In this respect, the phase change is, of course, 
very rapid, as it occurs during roughly a month. Towards the 
end of this segment, another phase jump is observed. The phase 
change is only A</> « 0.2 and based on the CPS analysis, it occurs 
gradually during ~ 20 days. 

In the beginning of SEG5 (HJD = 2451144.0394) two 
active longitudes reside at opposite sides of the star. Both the 
CF- and CPS -method results show that the regions are not com- 
pletely stable, but move in phase (upper right panel of Fig. [8]). 
This is also evident from the evolution of the light curve dur- 
ing this segment (Fig.©. At HJD ss 2451 170 there is a flip-flop 
switch between the spot regions and the former secondary min- 
imum becomes the new primary minimum. The CPS-analysis 
also reveals that at HJD « 2451280 one of the spot regions 
starts drifting towards the other one and seems to pass it around 
HJD a 2451300 (panel (g) in Fig. |5). The drift continues un- 
til HJD w 2451330. This could be explained by the fact that 



differential rotation will make surface features at different lat- 
itudes or with different anchor depths move with different an- 
gular velocity. This drift is not apparent in the CF-analysis, as 
expected, because the model is not supposed to register such fast 
changes. Instead we see a phase jump of A<p « 0.3 occurring 
after HJD « 24515300. From the plots of the individual light 
curves of the CPS-analysis we clearly see that the changes are 
actually gradual, and not just an effect of "interpolation" over an 
abrupt change (Fig. [5). Panels (h)-(j) of Fig. [5] show that there 
are clear connections between the light curve mean, amplitude 
and period. Especially interesting is the initial correlation be- 
tween A and M, which is evolved to a loop in the (A,M) dia- 
gramme. Correlations between neighbouring points in this kind 
of diagramme are expected, since the se are not independent mea- 
surements (see lLehtinen et alll201 ll) . But the loop indicates that 
there really is a more complex connection between A and M. 

The next detectable phase shift occurs in SEG8 (middle 
left panel of Fig. [8), during which the primary minimum sud- 
denly drops in strength, while a secondary minimum roughly 
A(p » 0.5 apart gains in magnitude, and during a few tens of 
days, the activity becomes concentrated to the location of the 
former secondary minimum. Although this event fulfils the cri- 
terion for a flip-flop, the new primary active longitude does not 
appear to become completely stable. Another minimum emerges 
at HJD « 2452380 nearby it, gains in magnitude, and finally the 
two minima seem to merge. 

In the beginning of SEG11 (middle right panel of Fig. [SJ, 
quite an abrupt phase shift, of nearly A0 = 0.5, can be seen. This 
event, however, is very close to the beginning of the data set and 
it is not detectable from the CPS analysis, probably due to the 
low light curve amplitude. 

In the beginning of SEG12 (lower left panel of Fig. [8), the 
primary minimum is located at — 0.6 roughly for the first 60 
days of observations. At HJD 2453770 a secondary minimum 
emerges at <p = 0.9, and an exchange of activity levels occurs 
between the active longitudes. At the same time, the former pri- 
mary minimum slowly drifts to <p — 0.2. This event takes over 
100 days 

In the beginning of segment SEG15 (lower right panel of 
Fig. [8]), both the CF and CPS-analysis reveal two active regions 
separated by ~ 0.5 in phase. Shifts in the strengths of these seem 
to occur at HJD « 2454830 and HJD * 2454870. The two ac- 
tive longitudes seem to merge and form a common minimum at 
HJD * 2454940. 

In conclusion, the CF and CPS analysis disclose complex 
phase behavior in six of the analysed segments, namely SEG4, 
SEG5, SEG8, SEG1 1, SEG12, and SEG15. The shifts in the pri- 
mary minima in these segments can be explained by: 

- Flip-flops as defined earlier in this section (SEG4, SEG5, 
SEG8, SEG1 1, SEG15; denoted "ff" in Table© 

- Phase jumps of A<p < 0.4 (SEG4, SEG12; denoted "phj") 

- Two drifting active regions (SEG5, SEG12; denoted "dr") 

- Merging of two active regions (SEG8, SEG15; denoted 
"mr"). 

The merging of the active regions do not necessarily mean a 
physical merger, but rath er that the two regions form a common 
minimum. As showed bv lLehtinen et alj d201 ll) . there is a mini- 
mum phase difference for spot regions under which these cannot 
be observed to cause separate minima. 

The five flip-flops occur at times (in years) t ss 1998.2, 
1999.0, 2002.0, 2005.0 and 2009.1. Thus, the interval between 
these events is 0. 8^4. 1 years and no clear periodicity can be seen. 
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Fig. 8. The photometric minima from the CF and CPS analysis of segments from left to right SEG4, SEG5, SEG8, SEG1 1, SEG12 
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Furthermore, the division between flip-flops and phase jumps is 
not completely clear, since there is a "grey zone" of A<p * 0.4 
phase shifts. We cannot rule out that the flip-flops and phase 
jumps play a role in the possible activity cycle, but the present 
data is not sufficient for any definite conclusions. 



4.3. Differential rotation 

In all segments we see drifts of the active regions. These indi- 
cate that the spots are not rotating with constant angular velocity. 
From the independent CPS period estimates we got the weighted 
mean period P w + AP W a 2<?3975 + $0123. If the variations in 
the photometric period were caused by differential rotation, we 



could estimate the differential rotation coefficient with the pa- 
rameter Z = (|Jetsulll993h . For FK Com, we got the value 
Z « 0.0308, which would correspond to a differential rotation of 
ACl > 4.6°/d. 

An alternative way to estimate the differential rotation is to 
study the drifts of active regions. During segment SEG5 the 
separation of the two active regions changed from -0.5 to 0.2 
in 60 days. This would correspond to a differential rotation of 
AQ « 4.2°/d between the two spot structures. This is, as it 
should be, less than the value derived from the variations in the 
rotation period, but st ill more than twice the value derived by 
Korhonen et al. (2007) assuming a solar differential rotation law 
(Eq.D. 
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Fig. 9. All reliable light curve fits for SEG5. 



4.4. Comparison with previous results 

Our analy s is is p artly based on the same data as used by 
Olah et alJ (|2006) and covers the same time as the D oppler 
images by Korhonen and collaborators (Korh onen et al.l l2007, 
2009a b). In Fig.[8]we already compared our results with simul- 
taneous Doppler imaging maps. In Fig. [10] we compare the pri- 
mary minima from the CF analysis with Doppler images from 
the years 1993-2008. This figure illustrates the rapid changes in 



the spot activity of FK Com. In general the consistency between 
the CF analysis and Doppler images is satisfactory. However, 
one must take into account that the photometric minimum is a re- 
sult of the integrated effect of several spot regions, including the 
visibility effect. In calculating the longitudinal Doppler imaging 
stripes, visibility and limb darkening was not taken into account. 

Interestinglv lOlah et al.1 (120061) also recognised the complex 
phase behaviour in SEG5. They conclude that a flip-flop oc- 
curred during this segment, while we detect both a flip-flop and 
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Fig. 10. Longitudinal slices of Doppler images and photometric minima from the CF analysis. 



Table 3. Summary of the flip-flop like events in FK Com: Time, 
duration and type of event. 



Segment 


HJD 


year 


At [d] 


Event type a 


SEG4 


2450870 


1998.2 


35 


ff 


SEG4 


2450960 


1998.4 


20 


Phj 


SEG5 


2451170 


1999.0 


10 


ff 


SEG5 


2451270 


1999.3 


60 


cli- 


SEG8 


2452270 


2002.0 


20 


ff 


SEG8 


2452420 


2002.4 


40 


mr 


SEG11 


2453360 


2005.0 


30 


ff 


SEG12 


2453790 


2006.1 


30 


Phj 


SEG12 


2453790 


2006.1 


110 


cli- 


SEG15 


2454870 


2009.1 


30 


ff 


SEG15 


2454940 


2009.3 


50 


mr 



1 ff = flip-flop, phj = phase jump, dr : 
merging spot regions 



drifting spot regions, mr : 



drifting spot regions. Furthermore, they also detect a phase jump 
in SEG3. In our analysis this would be the "bump" visible in Fig. 
[2] Similar small bumps can be seen in several of the segments, 
which is why we did not focus on this case. 



5. Discussion 

Our close study of six flip-flop like events clearly shows that 
we are not dealing with a singular phenomenon. The shifts of 
the active longitudes can be a result of both abrupt and gradual 
changes. In one case the apparent flip-flop can be explained by 
two spot regions moving with different angular velocity and even 



passing each other. The apparent flip-flop can thus be a conse- 
quence of differential rotation in two alternative ways: 

- The spot latitudes change and the spot or spot groups drift 
with respect to each other because of surface differential ro- 
tation. 

- The anchor depth of the spot or spot groups change and the 
differences in the angular velocity is caused by the depth de- 
pendent differential rotation. 

Furthermore, the situation will be complicated by rapid spot 
evolution. At times there are significant changes in the light 
curves within ~ 10 days. This time would probably be much 
shorter with light curves with denser spacing and higher accu- 
racy, e.g. satellite observations. The changes in the light curves 
of FK Com cannot be explained just by differential rotation of 
a steady spot model (see e.g. Fig. [3J. In reality we may thus be 
witnessing a combination of several effects: Rapid spot evolution 
combined with differential rotation in both depth and latitude, 
spiced with a possible dynamo wave. 

We find that the differential rotation in FK Com is at least 
about AQ m 4.6 deg /d. Which would correspond to a differen- 
tial rotation coefficient of k ss 0.03. This is roughly three times 
larger than the value suggested by iKorhonen et alJ d2007l) and 
about one sixth of the solar value. Our new value is in fair agree- 
ment with the observational a nd theoretical consensus (see e.g. 
iKitchatinov & Rtidigerl [l 999) of differential rotation diminish- 
ing proportional to 



AQ 



(8) 



where n * 0.8 - 0.9. The rotation period of the Sun being 
roughly an order of magnitude slower than that of FK Com, 
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indicates roughly 7 times weaker differential rotation for FK 
Com. However, we cannot conclude that this is a measure of 
the surface differential rotation, as it could also reflect geo- 
metri c properties of the l arge-s cale dynamo field, as described 
by iKorhonen & Elstneri d201 lb . and i nclude a signal from a 
possible azimuthal d ynamo wave (e.g. Srause & Radler 1980; 
ILindborg et aT]l201 ll) . 

The finding of an active longitude period, which is slightly 
longer than the mean rotation period indicates the presence of 
an azimuthal dynamo wave. In the case of FK Com this would 
rotate s lower than the star itself, as opposed to the RSCVn star 
II Peg (ILindborg et alj[2QTlt iHackman et al.ll20"TTl \m% . Even 
though a clear long-term active longitude period can be detected, 
the short term spot evolution seems fast and random. Thus, it is 
not as easy to follow active longitudes over gaps in the data as 
in the case of e.g. II Peg. In this se nse FK Com resem bles some 
other single stars, e.g. HP 1 16956 dLehtinen et al.ll201lh and LQ 
Hya dLehtinen et al.l 20l2h . In general, the active longitudes of 
close binary stars may be more regular due to the tidal effects. 

We also note that there is a clear tendency of having two ac- 
tive structures with a phase difference of considerably less than 
0.5. This tendency can also be seen in the distribution of the min- 
imum phases folded by P^: The maximum in Fig. [7] is divided 
into two peaks with a separation of roughl y 0.2. The same was 
already noted in earlier studies. E.g. Korho nen et al.l d2009bl) re- 
ported a pair of spots in 2008 with a phase difference of 0.25. 
This could resemble the pattern of leading and trailing spot pairs 
seen in much smaller scale on the Sun. 

Concerning the possible existence of an activity cycle our 
study does not bring much conclusive evidence. The mean mag- 
nitude shows variability which is compat ible with the cycle 
of 4.5-6.1 years found by lOlah et aD d2009l) . Naturally, when 
two active longitudes coexist, the amplitude tends to be lower. 
However, we do not see any other connection between the 
switches in the active longitudes and the long-term photometric 
variability. Neither can we detect any regularity in the flip-flop 
like events of FK Com. 
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